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Abstracts 

This paper proposes a conductance probe wtth 
vertical multiple electrode array (VMEA) which can be 
used in ofl-iu- water two-phase pipe flow to measure the 
water volume fraction aad axial velocity* The electrodes 
are flush mounted axiafly on the inside wall of an 
insulating duct Using finite element analysts method, the 
evaluating Indicators of uniform degree^ spatial sensitivity 
and Information volume were used to Investigate the 
performance of the VMEA. The optimized VMEA has 
enhanced spatial sensitivity and capability of obtaining the 
distributed Information from two-phase flow. Doe to the 
significant difference of conductance between oil and 
water, the VMEA can also be used to measure gas-liquid 
and liquid-solid two-phase How with conductive 
continuous phase. The evaluating indicators used In this 
paper can be used extensively to Investigate other types of 
conductive probe. 
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1. Introduction 

Oil-water two-phase flow is a common situation in 
oil logging and chemical engineering, so die 
measurement of oil-water is significant to the industrial 
process. A frequently used method of the measurement 
of two-phase flow is to set-up a theoretic model 
between volume fraction and the two-phase mixture 
impedance and apply cross-correlation algorithm to 
obtain the axial velocity of the mixture. The major 
advantages of this method are no interference to the 
flow, high frequency response, satisfying accuracy and 
the simplicity of industrial application. A 
comprehensive review of impedance method applied to 
two-phase flow was presented by Ceccio et ai. [ \ 

With reference to the impedance method, a plate 
electrodes probe was employed by Coney ^ to 
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investigate the theory behavior of flat electrodes for the 
measurement of liquid film thickness in two-phase flow. 
Dykesteen et ai. (3 ' described a theoretic model used to 
relate the measured impedance and void fraction in 
three-component flow. Kytmdaa et al w applied the 
guarded arc probe to validate the consistency of 
cross-correlation velocity and kinematic wave velocity. 
Tsochatzidis et al. 151 developed die mathematical model 
of Coney M to ring electrode system. Liu Xingbm 161 
presented a four-rmg-electrode conductive sensor for 
the measurement of water volume fraction. Chul-Hwa 
Song et al. 171 studied die void fraction using guarded arc 
probes and eliminated the error by introducing a 
reference probe. Fossa 1,1 investigated the performance 
of plate electrode and ring electrode. Devia&Fossa 191 
designed the shape of probe along pipe circumference 
to obtain the linear dependence between the void 
fraction and the measured mixture conductance. 

As far as axial velocity of two-phase flow is 
concerned, LucasA Walton' 103 applied a kinematic wave 
theory to measure the vertical upward bubbly two-phase 
flow. Lucas et al. [lu 121 embedded six ring electrodes on 
the surface of a mini insulating probe to measure 
liquid-solid two-phase flow and validate the ERT 
system. The ring conductance electrodes have been used 
by Lueas&Jin 1 ' M| to measure oil-in-water two-phase 
flow using the modified kenematic wave theory. An 
axial six-electrode sensor was used by Hu Jiohai et al. 
t,$I to obtain the water volume fraction and axial 
velocity of two-phase flow simultaneously 

Currently the plate electrode is applied to 
rectangular pipe and the arc or ring electrode is 
frequently used in cylinder pipes. The ring electrodes 
are mounted along the circumference of the pipe which 
has better spatial sensitivity. The performance of ring 
electrode array is decided by the electrode geometry, 
including the electrode height and the distance between 
the electrodes. Simultaneous measurement of volume 
fraction and velocity has been achieved using ring 
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six-electrode array but the measurement of void fraction 
was attached to the cross-correlation velocity 
measurement electrodes. It is difficult to confirm its 
significance of optimization. _ The geometry of 
conductance probes used now is based on experience or 
simple analysis. An important aspect, the electric field 
edge effect of rings, is always not considered which 
could lead to unexpectable error. Therefore, there exist 
the problems of geometric optimization and electrode 
arrangement. There is still a lack of comprehensive 
study on the theoretic evaluating indicators, so it is 
necessary to improve the existing VMEA. 

In order to tackle the problems discussed above, 
mis paper proposes a probe with new geometry which 
can be employed to measure oil -in- water two-phase 
pipe flow. The optimum geometry figuration has been 
obtained by investigating the uniform degree, spatial 
sensitivity and information volume of the VMEA in a 
1 25mm diameter cylinder pipe. 




(a) (b) 

Figure 1 Geometry of the VMEA probe 



2. Geometry of the VMEA for oil-in-water 
two-phase flow measurement 

As shown by figure 1(a), the new VMEA consists 
of eight conductance rings which are flush mounted on 
the inside wall of the pipe. Ej and E 2 are exciting 
electrodes. C r C 2 and C 3 -C 4 axe two up and down 
stream sensors denoted A and B respectively. H r H 2 are 
volume fraction sensor denoted C. All representing 
parameters used in this paper are shown in figure 1(b). 

As electrodes E t and Eq are excited by electrical 
source, the fluctuating signals across measuring 
electrodes (C,-C 2 , CrC 4 and H r H 2 ) reflect the 
conductance distribution of the two-phase flow. The 
signal from sensor C contains a mass of distributed 
information of oil-b-water two-phase flow, from which 
the water volume fraction can be obtained by a built 



model of water volume fraction. When two-phase 
mixture flows steadily inside the pipe, the signal from 
sensor A is similar to that from sensor B, on the basis of 
which the axial velocity can be obtained using 
cross-correlation algorithm. In this study, oil-m-water 
dispersed bubbly* flow is assumed to find the optimum 
configuration of the VMEA. 

3. Optimization of cidting electrodes 

3*1 The analysis of sensitive field of the VMEA 

The measuring principle of the VMEA is due to the 
theory of sensitive field. In order to prevent the 
corrosion of those electrodes and the electrolysis of 
two-phase mixture in the vicinity of electrodes, the 
titanium alky material and alternating voltage are 
applied. It is usually decided to keep the exciting 




Figure 2 The meshed 2D and 3D finite element model 



frequency of the applied voltage in the range of less 
than 100 KHz to eliminate the role of capacitance. For 
the inner pipe diameter is much small in comparison to 
the wavelength of the sensitive field, the electric field 
can be modeled as time-invariant, so Laplace equation 
Z*U=0 of static field can be used. 1 

The 2D and 3D models are meshed with triangle 
and tetrahedron respectively, which are shown as figure 
2. Material properties of electrodes and fluid mixture 
are taken into account The both ends of pipe should be 
assumed as infinite reference plane to eliminate the 
calculation error. One point worth noticing from the FE 
meshing is that the mesh cells should be fine enough to 
figure out the interracial effect between electrodes and 
fluids. 

The simulation is performed for all electrodes 
thickness of 1mm, exciting electrodes El and E2 
positioned at z=±2t/2 and ±10V voltage supplied The 
electric field distribution is shown as figure 3. As can be 
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seen, in the region adjacent to E, ami E 2l the distribution 
of vohage U is very peaky and approximately equal to 
the applied voltage; while in the range far from E, and 
E 2 , the U is tend to be linear. Figures 3(a) and 3(b) show 
that in the vicinity of £| and E 2 , both of the current 
density J f and J t vary sharply. Within the scope between 
exciting electrodes, the J r and J t distribute uniformly 
and tend to be constant Here this region is defined as 
uniform field region. Previous works have presented 
measuring electrodes should be set in this region, but 
further investigation in this study show that the 
optimum probe configurations are far beyond the 
previous works. 

3.2 The effect of H* on sensitive field 

For uniformly dispersed bubbly oi Lin-water flow, 
the sensitive field distribution is mainly effected by the 
height and distance of exciting electrode. 

The axial length of the uniform field section will 
be used following, of which quantization error is 
confined in the range Less than 1%. The current density 
near the pipe wall and axis is denoted J w and J A 
respectively. In order to investigate the uniform degree 
(UD) of the sensitive field along the radial and axial 
direction, the radial UD and axial UD are defined and 
denoted X^ and JW Xud and Yud can be calculated 
using the equations below: 

G/ r rH/J)/2 Z 4 

where Luf is the length of the uniform field region. 
The curves of J and Xyo are shown as figure 4 ( Very 
little variation of Yud is observed, so the Yud is not 
shown here). It can be seen that J increases with ft but 
Xvd decreases. As known, the polarization will 
increases with the expansion of interface between 
electrode and mixture, therefore 5mm is best 

3 J The effect of Z, on sensitive field 

Several evaluating indicators of sensitive field are 
calculated and shown as figure 5. As can be seen, the 
uniformity of field distribution can be improved with 
the increase of Z n but the current density J declines 
simultaneously. It is difficult to determine the value of 
Z* so further works are performed following. 

Since the edge effect of electric field, some part of 
the potential distribution is outside of the region 
between electrode E ( and electrode E 2 . As a result, the 
output of the VMEA will unavoidably include the 



ineffective distributed information of the fluid body. So 
the further work must be carried out to minimize the 
ineffective information. 

4. Optimization of volume fraction electrodes 

4.1 Spatial sensitivity distribution of volume 
fraction electrodes 

The water volume fraction can be obtained by the 
equivalent mixture conductance which contains the 
fluid information. The capability of getting information 
from fluid region is decided by the spatial sensitivity of 
measuring electrode, so the investigation of spatial 
sensitivity and the capability of obtaining mformarion is 
necessary. 

Lucas et a/. 1 11,151 investigated the spatial sensitivity 
of die axial six-electrode local probe which was inserted 
into the pipeline. The Lucas spatial sensitivity is used 
here. It is assumed that a small nonconducting fluid 
particle is positioned in the pipe with its center at 
coordinates (r,z) in r-z plane. 6U(r,z) is the change in 
voltage drop on the measuring electrode pair. Then the 
sensitivity of the measuring sensor ¥{r t z) is defined as: 



, . 6U(r % z) 



(2) 



where (SU)^ is the maximum SU(r t z\ then the 
sensitivity of sensor C can be calculated and denoted 
Vfcz). The calculated distribution of Yd**) is shown 
m figure 10. It can be seen that the Vdfrz) is high 
relative to that round the pipe axis in the region close to 
electrodes H, and H 2 The plot shows four distinct peaks 
alongside of the pipe wall. Clearly the spatial sensitivity 
is symmetrically distributed. One of the outboard peaks 
is between the electrodes H, and Ei whilst the other is 
between electrodes H 2 and E*. Both inboard peaks are 
between electrodes H, and H 2 and close to the 
corresponding electrode. 

In order to explicitly analyze the effect of the 
probe configuration on die spatial sensitivity 
distribution, here define the sensitivity along r axial as 
eigen sensitivity function denoted vdKz)^ and 
calculated for a wide range of electrode geometries. As 
figure 6 shown, ycfcQ) ^ determined by Z, and Z*. The 
plot exhibits significant difference between ycfoO) 
close to the pipe wall and round the pipe axis as Z t is 
equal to 100mm. The approximately horizontal curve 
appears when 3,=200mm and Z, -300mm, which means 
the electrode pair C can obtain approximately equal 
information from various fluid region. 
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42 Information of volume fraction electrodes 

The capability of obtaining fluid information 
should be enhanced when the spatial sensitivity 
distributed uniformly. Primarily, the fluid distributed 
information should be quanutated The information 
volume is defined as: 

V l (r t z)=j V rds=l^(r 3 z)drdz (3) 

3 S 

where V, ts the obtained total information volume 
of measuring electrode and where S is the region of 
measuring electrode pair. It can be seen from figure 10 
that ¥cfor) covers some part of the region outside the 
sensor C, with the result that some ineffective 
information is included unexpectively. The investigation 
of effective information should be performed. Effective 
information can be expressed as: 

y a (r 9 t) = 2 J* J*^ y,(r, z)drdz (4) 

where V& is (he effective information volume and 
the integral region is the field between the measuring 
electrodes. The concentration of V a in K ; is calculated 
using the expressions as below: 

C„(r,z) = ^^xI00 (5) 

The information volume curves of sensor C are 
shown as figure 7, where V J C> V D _ c *nd C DJ q denote the 
total information volume, effective information volume 
and concentration of effective information respectively. 
It can be seen that K / C and V £JC has always a maximum 
value and C a _ c increases when they arc «»i^iiaf^ for a 
wide range of electrode geometries. For measuring 
electrodes, the criterion to determine the value of 2* 
should be that it can get as most effective information 
as possible from the measured fluid region. 

It is worthwhile to notice that there is little 
interaction between volume traction sensor C and cross 
correlation sensors A and B, so the investigation of 
sensors A and B can be performed without modeling 
volume fraction sensor. 

5. Optimization of crost-co relation electrodes 

Lucas et al) 12 * described that the 'sensing volume* 
should be compact enough to keep high accuracy of the 
transition spacing of the fluid body flowing from 
upstream sensor A to downstream sensor B which is 
denoted Zci^see figure 1(b)). In order to get accurate 
transit time, the C a of sensor A and sensor B should be 
the maximum value. 



According to the spatial sensitivity defined in 
equation (2), the spatial sensitivity of sensor A and 
sensor B arc denoted !%2) and Yfaz) respectively 
and calculated in a similar way to Figure] 1 

exhibits the symmetrical relation between ^fcz) and 
VXnz). The peak of ¥Ar>2) is approximately midway 
between electrodes C, and C 2 whilst the peak of yj(r P r) 
is approximately midway between electrodes C 3 and CV 
The peak value is defined as the eigenvalue to 
investigate its variational regulation. The calculated 
eigenvalue distribution is shown as figure 8 f of which 
the influencing factors are Zee, Z t and Dec* It can be 
seen that as the position of sensors A and B is towards 
the center region of the VMEA, V a (kz) and tFfcfrz) 
decreases simultaneously, hence sensors A and B should 
be not too close. Since the position of V H =100(or 
t^IOO) always appears as Dec is equal to 20mm, an 
important conclusion can be drawn apparently from 
figure 8 is that the optimum value of Dec should be 
20mm. 

As figure 1 1 shown! some part of f A (nz) is outside 
of the effective region of sensor A whilst some part of 
Vfaz) is outside of the effective region of sensor B, 
which demonstrates that ineffective mformatton exists 
in the total information of the output cross sensors A 
and B. In order to get accurate transit time, the 
concentration of effective information of sensor A and 
sensor B is investigated for D«r20mm{Figure 9). As 
the position of sensor A (or B) is towards electrode E, 
(or E*) and far from the position z=0, Q,^ and C a B 
decrease, where C^and C^are the concentration of 
the effective information of cross-correlation electrodes 
A and B. As Z, increases, Csu and decrease as 
well All the evaluating indicators should be taken into 
account as the optimum probe configuration is to be 
determined. 

As demonstrated by figure 7, there always exists a 
maximum of Vu_ c for each group of Z k and Z 4 and the 
difference between the every two maximums is distinct 
Figure 6 shows the approximately horizontal IrHnO) 
curve appears when Z, is over 200mm. The distinct 
difference between the sensitivity of sensor A and 
sensor B can also be seen in figure8. Taking above 
discussions into account, 250mm should be the best 
value of Z^ then Z A =170mm can be obtained from 
figure 7. Figure 10 also shows 55mm is the preferable 
Caji or Oca, the corresponding Zee is 110mm, The 
spatial sensitivity distribution of the optimized VMEA 
is shown in figure 10 and figure LI. 
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6. Conclusions 

Based on Laplace problem of electric field of the 
VMEA and the systematic evaluating indicators, the 
geometry of vertical multiple electrode conductive array 
is optimized. The characteristics of the optimized 
VMEA can be concluded as below: 

The volume traction sensor is separated from 
cross-correlation sensors, which can enhance the spatial 
sensitivity and capability of acquiring equal information 
from fluid regions. 

The cross-correlation sensors have maximum 
concentration of effective information and small 
measuring space, which guarantee that accurate 
measurement of fluid velocity can be achieved. 

The optimized VMEA can measure volume 
fraction and the fluid velocity of oiMn- water two-phase 
flow simultaneously. 
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D = 125mm, Z e = 90mm, H e = 5mm, H m = 3mm, a w = 0.0 ISm' 1 , a n = 2.3 x 10 7 Sm\ a = 1 

cr^, - water conductivity, c ri - conductivity of titanium, a- water volume fraction 
U - voltage distributed in the VMEA, Jr - the component of current density in radial direction 
Jz - the component of current density in axial direction 
Figure 3 Spatial sensitivity distribution of the volume fraction sensor C 
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Figure 8 Eigenvalue distributions of cross-correlation sensors 
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